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""Fully  polarimetric  scattering  properties  of  earth  terrain  media  are  studied  with 
three-layer  random  medium  model.  The  top  scattering  layer  is  modeled  as  an  isotropic 
random  medium  which  is  characterized  by  a  scalar  permittivity.  The  middle  scattering 
layer  is  modeled  as  an  anisotropic  random  medium  with  a  symmetric  permittivity  ten¬ 
sor  whose  optic  axis  cam  be  tilted  depending  on  the  preferred  alignment  of  the  embedded 
scatterers.  The  bottom  layer  is  considered  as  a  homogeneous  half-space.  Volume  scatter¬ 
ing  effects  of  both  random  media  are  described  by  three-dimensional  correlation  functions 
with  variances  and  correlation  lengths  corresponding  to  the  strengths  of  the  permittivity 
fluctuations  and  the  physical  sizes  of  the  inhomogeneities,  respectively.  The  strong  fluc¬ 
tuation  theory  is  used  to  derive  the  mean  fields  in  the  random  media  under  the  bilocal 
approximation  with  similarities  of  the  dyadic  Green’s  functions  properly  taken  into  consid¬ 
eration.  With  the  discrete  scatterer  concept,  effective  permittivities  of  the  random  media 
are  calculated  by  two-phase  mixing  formulas.  Then,  the  distorted  Bora  approximation  is 
used  to  calculate  the  covariance  matrix  which  describes  the  fully  polarimetric  scattering 

properties  of  the  remotely  sensed  media.  The  polarimetric  information  is  useful  in  the 
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To  observe  polarimetric  scattering  directly  from  earth  terrain  media  such  as  veg¬ 
etation,  meadow,  and  ice,  the  three-layer  configuration  is  first  reduced  to  two-layer  by 
removing  the  top  layer.  Such  media  exhibit  reciprocity  as  experiment  ally  manifest  in  the 
close  proximity  of  the  measured  backscattering  radar  cross  sections  ovh  and  c»v  and  the¬ 
oretically  established  in  the  random  medium  model  with  a  symmetric  permittivity  tensor. 
In  this  case,  the  covariance  matrix  is  a  3  by  3  matrix  with  the  lower  triangle  containing  el¬ 
ements  which  are  complex  conjugates  of  the  corresponding  elements  in  the  upper  triangle. 
The  results  obtained  with  the  two-layer  configuration  show  that  polarimetric  information 
carried  by  the  complex  correlation  coefficient  (p)  is  helpful  in  identifying  isotropic  versus 
anisotropic  media  due  to  the  distinctive  characteristics  that  each  of  the  two  media  imposes 
on  p.  For  an  isotropic  medium,  the  value  of  p  is  close  to  1  over  the  range  of  incident  angles 
under  the  single  scattering  approximation.  Whereas  for  an  anisotropic  medium,  the  mag¬ 
nitude  of  p  attains  a  maximum  value  at  a  particular  incident  angle  6P  and  decays  as  the 
incident  angle  departs  from  9P.  Furthermore,  the  correlation  coefficient  conveys  informa¬ 
tion  about  the  optic-axis  tilted  angle  of  the  anisotropic  random  medium.  When  the  optic 
axis  is  not  tilted,  9P  is  at  normal  incident.  As  the  optic  axis  becomes  tilted,  9„  deviates 
away  from  normal  incident.  It  should  also  be  noted  that  the  tilt  of  the  optic  axis  is  di¬ 
rectly  related  to  the  non-zero  depolarization  terms  in  the  covariance  matrix  which  becomes 
very  small  for  an  untilted  anisotropic  or  isotropic  medium.  Thus,  the  fully  polarimetric 
information  is  important  in  the  remote  sensing  of  earth  terrain. 

The  three-layer  random  medium  has  the  capability  of  accounting  for  polarimetric 
scattering  from  earth  terrain  media  under  the  effects  of  weather,  seasonal  variation,  and 
atmospheric  conditions  such  as  forest  under  mist,  meadow  under  fog,  and  ice  under  snow. 
The  effects  on  polarimetric  wave  scattering  due  to  the  top  layer  are  identified  by  comparing 
the  three-layer  model  results  with  those  obtained  from  the  two-layer  model.  Consider  first 
the  case  of  a  low-loss  scattering  layer  covering  the  tilted  anisotropic  scattering  layer.  The 
enhancement  of  the  radar  returns  due  to  dry-snow  cover  on  top  of  first-year  sea  ice  observed 
in  the  experimental  data  can  be  explained  using  the  three-layer  configuration.  Also,  the 
low-loss  top  layer  can  give  rise  to  the  oscillation  on  the  radar  returns  as  a  function  of 
incident  angle  due  to  the  boundary  effect.  The  oscillation  can  also  be  seen  on  the  real  and 
imaginary  parts  of  the  correlation  coefficient  (p).  Interestingly,  the  magnitude  of  p  does  not 
show  the  oscillation  while  clearly  retaining  the  similar  characteristics  as  observed  directly 


from  the  two-layer  configuration.  Thus,  the  correlation  coefficient  can  carry  information 
from  both  the  covering  low-loss  isotropic  layer  and  the  lower  tilted  anisotropic  layer  in  a 
rather  distinctive  manner.  As  the  thickness  of  the  top  layer  increases,  the  value  of  p  at 
a  given  incident  angle  increases  meaning  that  the  top  layer  can  weaken  the  anisotropic 
effect  of  the  lower  layer.  This  masking  effect  is  clearly  seen  when  the  top  layer  becomes 
lossy.  In  this  case,  the  radar  returns  are  diminished  instead  of  being  enhanced.  In  this 
application  of  the  random  medium  model  to  polarimetric  remote  sensing  of  earth  terrain, 
the  encountered  media  are  reciprocal  and  can  be  characterized  by  symmetric  permittivity 
tensors. 

A  radar  clutter  model  is  used  to  simulate  fully  polarimetric  returns  for  a  stepped 
frequency  radar.  The  purpose  is  to  create  synthetic  site  dependent  clutter  signatures 
that  can  be  utilized  in  a  hardware-in-the-loop  test  system.  The  fully  polarimetric,  multi- 
frequency,  multi-incident  angle  random  medium  model  is  employed  to  generate  normalized 
backscatter  coefficients  of  terrain  clutter.  This  model  is  used  to  generate  the  polarimetric 
terrain  clutter  covariance  matrices  for  each  of  N  high  resolution  range  bins,  at  each  of  the 
M  discrete  frequencies.  The  random  medium  model  allows  us  to  include  the  effect  of  the 
terrain  local  incident  angle  on  the  clutter  covariance  matrix  elements.  In  the  simulation, 
we  assume  that  there  is  a  single  clutter  class  within  each  of  the  N  range  bins,  although  the 
depression  angle  may  vary  from  bin  to  bin.  The  covariance  matrices  are  decomposed  and 
multiplied  by  complex  Gaussian  noise  in  order  to  generate  the  normalized  electric  fields 
in  the  backscattering  direction  for  each  of  the  N  range  bins,  at  each  of  the  AT  discrete 
frequencies.  These  fields  are  then  coherently  added,  taking  into  account  the  effects  of  both 
terrain  elevation  and  range.  This  yields  a  single  frequency  polarimetric  return  that  a  radar 
would  measure  from  the  specified  terrain.  The  radar  return  for  each  of  the  other  discrete 
frequencies  is  calculated  in  a  similar  manner.  The  result  is  the  clutter’s  low  resolution  range 
polarimetric  profile,  i.e.,  the  backscattered  signal  response  within  the  beam  footprint  of 
the  radar  antenna.  Each  discrete  frequency  is  simulated  and  the  effects  of  shadowing  and 
overlay  are  taken  into  account.  The  simulation  produces  coherent  phase-history  clutter 
returns  which  can  be  coherently  superimposed  on  the  target  phase-history  returns.  The 
combined  (or  clutter  only)  returns  are  processed  to  obtain  either  (1)  the  coherent,  high 
resolution  range  profile  or  (2)  the  noncoherent,  autocorrelation  range  profile. 
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Supervised  and  unsupervised  classification  procedures  are  developed  and  applied  to 
synthetic  aperture  radar  (SAR)  polarimetric  images  in  order  to  identify  its  various  earth 
terrain  components.  For  the  supervised  classification  processing,  the  Bayes  technique  is 
utilized  to  classify  fully  polarimetric  and  normalized  polarimetric  SAR  data.  Simpler 
polarimetric  discriminates,  such  as  the  unnormalized  and  normalized  magnitude  response 
of  the  individual  receiver  channel  returns,  in  addition  to  the  phase  difference  between  the 
receiver  channels  are  also  considered.  Covariance  matrices  are  computed  for  each  terrain 
class  from  selected  portions  within  the  image  where  ground  truth  is  available,  under  the 
assumption  that  the  polarimetric  data  has  a  multivariate  Gaussian  distribution.  These 
matrices  are  used  to  train  the  optimal  classifier,  which  in  turn  is  used  to  classify  the  entire 
image.  In  this  case,  classification  is  based  on  determining  the  distances  between  the  training 
classes  and  the  observed  feature  vector,  then  assigning  the  feature  vector  to  belong  to  that 
training  class  for  which  the  distance  was  minimum.  Another  processing  algorithm  based 
on  comparing  general  properties  of  the  Stokes  parameters  of  the  scattered  wave  to  that 
of  simple  scattering  models  is  also  discussed.  This  algorithm,  which  is  an  unsupervised 
technique,  classifies  terrain  elements  based  on  the  relationship  between  the  orientation 
angle  and  handedness,  or  ellipticity,  of  the  transmitted  and  received  polarization  state. 
These  classification  procedures  will  be  applied  to  San  Francisco  Bay  and  Traverse  City  SAR 
imagery,  supplied  by  the  Jet  Propulsion  Laboratory.  It  is  shown  that  fully  polarimetric 
classification  yields  the  best  overall  performance.  Also,  in  some  selected  areas  where  the 
observed  amplitudes  of  the  returns  are  quite  different  than  that  of  the  training  data, 
classification  techniques  not  based  on  the  absolute  amplitudes  of  the  returns,  e.g.,  the 
normalized  polarimetric  classifier,  produced  a  more  consistent  result  with  respect  to  the 
ground  truth  data. 

The  normalized  polarimetric  classifier  is  proposed  such  that  only  the  relative  magni¬ 
tudes  and  phases  of  the  polarimetric  data  will  be  utilized  to  discriminate  terrain  elements. 
For  polarimetric  data  with  arbitrary  probability  density  function  (PDF),  the  distance 
measures  of  the  normalized  polarimetric  classifier  based  on  a  general  class  of  normaliza¬ 
tion  functions  are  shown  to  be  equivalent  to  one  another.  The  normalized  polarimetric 
classifier  thus  derived  will  be  optimal  among  all  normalization  schemes,  when  the  system 
absolute  calibration  factors  are  common  to  all  polarimetric  channels.  Further  assuming  a 
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multivariate  complex  Gaussian  distribution  for  the  un-normalized  data,  the  distance  mea¬ 
sure  of  the  normalized  polarimetric  classifier  is  given  explicitly  and  turns  out  to  be  also 
independent  of  the  number  of  scatterers.  The  usefulness  of  the  normalized  polarimetric 
classifier  is  demonstrated  by  the  classification  of  trees  and  grass  in  the  experimental  data 
obtained  from  Lincoln  Laboratory.  The  classification  error  is  shown  to  be  the  smallest 
among  those  of  magnitude  ratio  or  phase  difference  classifications. 

A  three-layer  random  medium  model  is  developed  to  study  the  fully  polarimetric 
scattering  properties  of  earth  terrain.  The  top  layer  is  modeled  as  an  isotropic  random 
medium,  the  middle  layer  as  an  anisotropic  random  medium,  and  the  bottom  layer  as  a 
homogeneous  half-space.  Volume  scattering  effects  of  both  random  media  are  characterized 
by  correlation  functions  in  which  variances  and  correlation  lengths  describe  strengths  of 
permittivity  fluctuations  and  physical  sizes  of  embedded  inhomogeneities,  respectively. 
The  anisotropic  effect  of  the  middle  layer  is  attributed  to  specific  structure  and  alignment 
of  the  scatterers.  With  the  strong  fluctuation  theory,  the  mean  fields  in  the  random 
media  are  derived  under  the  bilocal  approximation  with  singularities  of  the  dyadic  Green’s 
functions  properly  taken  into  consideration.  With  the  discrete  scatterer  concept,  effective 
permittivities  of  the  random  media  are  calculated  by  two-phase  mixing  formulas.  Then,  the 
distorted  Born  approximation  is  used  to  calculate  the  covariance  matrix  which  describes  the 
fully  polarimetric  scattering  properties  of  the  terrain  and  is  used  in  radar  image  simulation 
and  earth  terrain  identification  and  classification. 

A  two-layer  random  medium  model  has  been  successfully  applied  to  polarimetric 
remote  sensing  of  earth  terrain  such  as  vegetation,  meadow,  and  ice  layer.  The  results 
obtained  with  the  three-layer  configuration  have  the  capability  of  accounting  for  polari¬ 
metric  scattering  from  earth  terrain  under  the  effects  of  weather,  seasonal  variation,  and 
atmospheric  conditions  such  as  forest  under  mist,  meadow  under  fog,  and  ice  under  snow. 
The  effects  on  polarimetric  wave  scattering  due  to  the  top  layer  are  identified  by  comparing 
the  three-layer  model  results  with  those  obtained  from  the  two-layer  model.  The  enhance¬ 
ment  of  the  radar  returns  due  to  dry-snow  cover  on  top  of  first-year  sea  ice  observed  in 
the  experimental  data  can  be  explained  using  the  three-layer  random  medium  model.  The 
theoretical  results  are  illustrated  by  comparing  the  calculated  covariance  matrices  with  the 
polarimetric  measurement  data. 
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There  is  a  considerable  interest  in  determining  the  optimal  polarizations  that  achieve 
maximum  contrast  between  two  scattering  classes  in  radar  polarimetric  images  for  the  pur¬ 
pose  of  terrain  discrimination.  In  this  paper,  we  present  a  systematic  approach  for  obtain¬ 
ing  the  optimal  polarimetric  matched  filter  which  produces  maximum  contrast  between 
two  scattering  classes,  each  represented  by  its  respective  covariance  matrix. 

To  accomplish  this,  we  derive  a  linear  weighting  vector  that  maximizes  the  expected 
power  return  ratio,  i.e.,  the  contrast  ratio  between  the  two  scattering  classes.  The  maxi¬ 
mization  procedure  involves  solving  an  eigenvalue  problem  where  the  eigenvector  yielding 
this  maxima  will  correspond  to  the  optimal  polarimetric  matched  filter.  Then,  through 
use  of  polarization  synthesis,  it  is  demonstrated  that  when  this  weighting  vector  is  uti¬ 
lized  to  process  fully  polarimetric  radar  images,  the  maximum  contrast  between  the  two 
respective  classes  results.  The  sub-optimal  problem  of  a  fixed  transmitting  polarization  is 
also  considered.  In  this  case,  the  received  polarization  is  optimized  such  that  a  maxima  in 
the  contrast  ratio  is  obtained  under  this  constraint.  To  exhibit  the  physical  significance  of 
this  filter,  we  transform  it  into  its  associated  transmitting  and  receiving  polarizations,  in 
terms  of  their  horizontal  and  vertical  vector  components. 

This  technique  is  then  applied  to  radar  polarimetry  obtained  from  the  Jet  Propul¬ 
sion  Laboratory.  It  is  shown,  both  numerically  and  through  the  use  of  radar  imagery, 
that  maximum  image  contrast  can  be  realized  when  data  is  processed  with  the  optimal 
polarimetric  matched  filter. 

For  active  and  passive  microwave  remote  sensing  of  sea  ice,  the  two-layer  uniaxial 
random  medium  model  is  applied  to  study  the  volume  scattering  and  anisotropic  effects 
which  are  attributed  to  embedded  brine  inclusions.  In  the  model,  the  correlation  function 
is  the  direct  link  between  the  electrical  behavior  and  physical  properties  of  embedded  brine 
inclusions.  The  spatial  distribution,  size,  elongated  structure,  and  preferred  alignment  of 
embedded  brine  inclusions  can  be  described  by  variances,  correlation  lengths,  and  the  form 
of  the  correlation  function.  We  have  extracted  a  correlation  function  of  exponential  form 
from  the  photograph  of  a  horizontal  thin  section  prepared  from  a  sample  of  artificially 
grown  saline  ice  that  closely  resembled  Arctic  congelation  sea  ice.  It  is  found  that  the 
extracted  correlation  lengths  are  consistent  with  the  published  average  size  of  brine  pock¬ 
ets.  With  the  application  of  strong  fluctuation  theory  and  the  bilocal  approximation,  the 
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effective  permittivity  tensor  is  derived  in  the  low  frequency  limit  for  am  unbounded  uni¬ 
axial  random  medium  with  two-phase  mixtures.  Using  the  extracted  correlation  length, 
the  effective  permittivity  tensor  is  computed  as  a  function  of  fractional  volume  of  brine 
inclusions  and  compared  with  the  in  situ  dielectric  measurements  at  microwave  frequencies 
of  4.8  and  9.5  GHz.  (Electromagnetic  theory,  sea  ice,  two-layer  unisonal  random  medium 
model,  volume  scattering,  anisotropic,  brine  inclusions,  Arctic  congelation  sea  ice,  strong 
fluctuation  theory,  bilocal  approximation,  effective  permittivity  tensor,  two-phase  mix¬ 
ture,  fractional  volume.)  0659  Random  media  and  rough  surfaces,  0669  Scattering  and 
diffraction,  0689  Wave  propagation 

A  three-layer  random  medium  model  is  developed  to  investigate  effects  of  snow 
cover  on  the  sea  ice  signature  for  microwave  remote  sensing.  The  volume  scattering  and 
anisotropic  effects  due  to  embedded  inhomogeneities  in  the  snow-covered  sec  ice  are  studied 
with  the  wave  theory.  The  snow  layer  is  modelled  as  an  isotropic  random  medium  and 
the  sea  ice  as  an  anisotropic  random  medium.  Volume  scattering  effects  are  caused  by 
granular  ice  particles  and  water  contents  in  snow  and  by  brine  inclusions  and  air  bubbles 
in  sea  ice,  respectively.  The  anisotropic  effect  is  attributed  to  the  elongated  structure  of 
brine  inclusions  with  a  preferred  alignment  between  ice  platelets.  In  the  random  medium 
model,  the  essential  quantity  is  the  correlation  function  which  contains  important  physical 
parameters  such  as  variances  and  correlation  lengths  for  characterizing  the  strength  of  the 
permittivity  fluctuation,  the  physical  size,  and  the  geometrical  structure  of  scatterers.  We 
h*ve  extracted  correlation  function  from  digitized  image  of  saline  ice  sample  grown  in  the 
outdoor  tank  at  CRREL.  The  calculated  correlation  lengths  are  consistent  with  reported 
average  size  of  brine  pockets. 

The  strong  fluctuation  theory  is  applied  to  account  for  the  distinct  permittivity 
difference  between  air,  ice  particles,  water,  and  brine  inclusions.  First  of  all,  singularities  of 
the  dyadic  Green’s  functions  for  snow  and  sea  ice  are  properly  considered.  Then,  the  mean 
fields  in  both  media  are  derived  by  the  bilocal  approximation.  With  the  discrete  scatterer 
concept,  effective  permittivities  for  dry  snow  and  first-year  sea  ice  are  computed  by  two- 
phase  mixing  formulas.  The  effective  permittivities  of  sea  ice  calculated  as  a  function  of 
fractional  volume  of  brine  inclusions  are  found  to  be  consistent  with  the  published  results 
obtained  from  slab  ice  samples.  For  calculating  effective  permittivities  of  wet  snow  and 
multi-year  sea  ice,  three-phase  mixing  formulas  are  used  since  wet  snow  consists  of  air, 
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ice  grains,  and  water  and  multi-year  sea  ice  is  composed  of  pure  ice,  brine  inclusions,  and 
air  bubbles.  Tinally,  the  bistatic  scattering  and  backscattering  coefficients  are  computed 
with  the  distorted  Born  approximation.  The  results  are  illustrated  by  comparing  the 
backscattering  coefficients  as  a  function  of  incident  angle  with  experimental  data  from 
controlled  field  measurements. 

A  generalized  K-distribution  is  proposed  to  model  the  statistics  of  the  fully  po- 
larimetric  returns  from  the  terrain  cover.  In  the  past,  K-distribution  has  been  used  suc¬ 
cessfully  to  characterize  the  intensity  distribution  of  single  polarization  returns.  Recently, 
fully  polarimetric  synthetic  aperture  radar  (SAR)  data  with  polarizations  HH,  HV,  and 
VV  have  been  proven  useful  in  remote  sensing  of  earth  terrain  and  in  this  paper  we  gener¬ 
alize  the  K-distribution  to  model  fully  polarimetric  terrain  radar  clutter.  The  generalized 
K-distribution  is  a  better  description  of  the  statistics  of  the  SAR  polarimetric  data  than 
the  Gaussian  distribution  and  would  be  useful  in  identification  and  classification  of  terrain 
types  in  the  SAR  polarimetric  data. 

In  general,  all  polarization  returns  of  each  single  polarimetric  measurement  are 
correlated  with  different  variances.  Therefore,  we  assume  an  n-dimensional  anisotropic 
random  walk  model  where  the  coordinate  components  of  each  step  are  characterized  by 
a  covariance  matrix  and  the  number  of  steps  or  scatterers  is  of  negative  binomial  distri¬ 
bution  with  parameter  a.  The  anisotropy  in  the  model  refers  to  the  fact  that  the  the 
covariance  matrix  is  not  proportional  to  an  identity  matrix.  A  generalized  K-distribution 
for  polarimetric  data  is  derived  when  the  average  number  of  steps  approaches  to  infinity. 
The  K-distribution  is  also  generalized  to  the  non-zero  mean  case,  which  can  be  used  to 
model  the  statistics  of  transmitted  electromagnetic  wave  through  atmosphere.  It  is  found 
that  if  a  zero-mean  K-distributed  random  vector  is  normalized  by  its  Euclidean  norm,  the 
joint  probability  distribution  of  the  normalized  quantities  is  independent  of  the  parameter 
a  and  is  the  same  as  that  derived  from  a  zero-mean  Gaussian-distributed  random  vector. 
The  results  are  illustrated  by  analyzing  the  normalized  intensity  moments  of  the  polari¬ 
metric  SAR  images  provided  by  the  Jet  Prop’ilsion  Laboratory  and  comparing  them  with 
the  generalized  K-distributed  polarimetric  model. 
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The  volume  scattering  effects  of  snow-covered  sea  ice  are  studied  with  a  three-layer 
random  medium  model  for  microwave  remote  sensing.  The  strong  fluctuation  theory  and 
the  bilocal  approximation  are  applied  to  calculate  the  effective  permittivities  for  snow  and 
sea  ice.  The  wave  scattering  theory  in  conjunction  with  the  distorted  Born  approximation 
is  then  used  to  compute  bistatic  coefficients  and  backscattering  cross  sections.  Theoretical 
results  are  illustrated  by  matching  experimental  data  for  dry  snow-covered  thick  first-year 
sea  ice  at  Point  Barrow.  The  radar  backscattering  cross  sections  are  seen  to  increase  with 
snow  cover  for  snow-covered  sea  ice,  due  to  the  increased  scattering  effects  in  the  snow 
layer.  The  results  derived  can  also  be  applied  to  the  passive  remote  sensing  by  calculating 
the  emissivity  from  the  bistatic  scattering  coefficients. 

The  remote  sensing  of  sea  ice  is  studied  with  the  two-layer  random  medium  model 
where  a  correlation  function  is  used  to  characterize  the  randomly  fluctuating  part  of  the 
permittivity.  Due  to  the  shape  and  distribution  of  the  brine  inclusions,  the  sea  ice  is  gen¬ 
erally  anisotropic.  We  assume  the  sea  ice  to  be  uniaxial  and  extract  correlation  functions 
from  a  photograph  of  an  artificially  grown  ice  sheet  sample.  The  sizes  and  distributions  of 
the  brine  pockets  are  related  to  the  correlation  lengths  and  the  shapes  of  the  correlation 
functions.  The  effective  permittivity  of  sea  ice  is  calculated  with  the  strong  fluctuation 
theory  assuming  a  two-phase  mixture  consisting  of  pure  ice  and  brine  inclusions.  The 
calculated  effective  permittivities  are  found  to  be  consistent  with  the  published  results 
obtained  from  slab  ice  samples. 

A  systematic  approach  for  the  identification  of  terrain  media  is  developed  using 
the  optimum  polarimetric  classifier.  The  covariance  matrices  for  various  terrain  cover  are 
computed  from  theoretical  models  of  random  medium  by  evaluating  the  full  polarimetric 
scattering  matrix  elements.  The  optimal  classification  scheme  makes  use  of  a  quadratic 
distance  measure  and  is  currently  applied  to  classify  a  vegetation  canopy  consisting  of 
both  trees  and  grass.  Experimentally  measured  data  are  used  to  validate  the  classification 
scheme.  Theoretical  probability  of  classification  error  using  the  full  polarimetric  matrix  is 
compared  with  classification  based  on  single  features  which  include  the  phase  difference  be¬ 
tween  the  VV  and  HH  polarization  returns.  It  is  shown  that  the  full  polarimetric  results  are 
optimal  and  provide  better  classification  performance  than  single  feature  measurements. 
The  application  of  this  classification  scheme  to  sea  ice  will  be  explored. 


10 


The  Mueller  matrix  and  polarization  covariance  matrix  are  described  for  polarimet- 
ric  radar  systems.  The  clutter  is  modelled  by  a  layer  of  random  permittivity,  described  by 
a  three-dimensional  correlation  function,  with  variance,  and  horizontal  and  vertical  corre¬ 
lation  lengths.  This  model  is  applied,  using  the  wave  theory  with  Born  approximations 
carried  to  the  second  order,  to  find  the  backscattering  elements  of  the  polarimetric  ma¬ 
trices.  It  is  found  that  8  out  of  16  elements  of  the  Mueller  matrix  are  identically  zero, 
corresponding  to  a  covariance  matrix  with  four  zero  elements.  Theoretical  predictions  are 
matched  with  experimental  data  for  vegetation  fields. 

Geophysical  media  encountered  in  nature  are  generally  mixtures  of  dielectrically 
different  materials.  When  a  dielectric  mixture  is  considered  macroscopically,  i.e.,  the  scale 
of  interest  is  much  larger  than  the  correlation  length  of  the  inhomogeneities  of  the  ma¬ 
terial,  it  can  be  assigned  an  effective  permittivity  that  relates  the  average  macroscopic 
displacement  and  the  electric  field.  This  effective  permittivity  can  be  calculated  through 
a  quasistatic  analysis  which  restricts  the  results  of  this  dielectric  mixture  theory  to  cases 
where  the  size  of  the  inhomogeneities  is  much  smaller  than  the  wavelength.  The  limitation 
resulting  from  this  is  that  the  effective  permittivity,  which  can  have  complex  values,  only 
includes  the  absorption  losses  but  not  the  scattering  losses. 

We  have  studied  the  volume  scattering  effects  of  snow-covered  sea  ice  with  a  three- 
layer  random  medium  model  for  microwave  remote  sensing.  The  strong  fluctuation  theory 
and  the  bilocal  approximation  are  applied  to  calculate  the  effective  permittivities  for  snow 
and  sea  ice.  The  wave  scattering  theory  in  conjunction  with  the  distorted  Born  approx¬ 
imation  is  then  used  to  compute  bistatic  coefficients  and  backscattering  cross  sections. 
Theoretical  results  are  illustrated  by  matching  experimental  data  for  dry  snow-covered 
thick  first-year  sea  ice  at  Point  Barrow.  The  radar  backscattering  cross  sections  are  seen 
to  increase  with  snow  cover  for  snow-covered  sea  ice,  due  to  the  increased  scattering  effects 
in  the  snow  layer.  The  results  derived  can  also  be  applied  to  the  passive  remote  sensing 
by  calculating  the  emissivity  from  the  bistatic  scattering  coefficients. 
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The  scattering  of  electromagnetic  waves  from  a  randomly  perturbed  peridic  surface 
is  solved  using  the  Extended  Boundary  Condition  (EBC)  method.  The  scattering  from 
periodic  surface  is  solved  exactly  using  the  EBC  method  and  this  solution  is  used  in  the 
small  perturbation  method  to  solve  for  the  scattered  field  from  a  randomly  perturbed  peri¬ 
odic  surface.  The  random  perturbation  is  modeled  as  a  Gaussian  random  process  and  the 
surface  currents  and  the  scattered  fields  are  expanded  and  solved  up  to  the  second  order. 
The  theoretical  results  are  illustrated  by  calculating  the  bistatic  and  backscattering  coef¬ 
ficients.  It  is  shown  that  as  the  correlation  length  of  the  random  roughness  increases,  the 
bistatic  scattering  pattern  of  the  scattered  fields  show  several  beams  associated  with  each 
Bragg  diffraction  direction  of  the  periodic  surface.  When  the  correlation  length  becomes 
smaller,  then  the  shape  of  the  beams  become  broader.  The  results  obtained  using  the  EBC 
method  is  also  compared  with  the  results  obtained  using  the  Kirchhoff  approximation.  It 
is  shown  that  the  Kirchhoff  approximation  results  show  quite  a  good  agreement  with  EBC 
method  results  for  the  VV  and  HH  polarized  backscattering  coefficients  for  small  angles 
of  incidences.  However,  the  Kirchhoff  approximation  does  not  give  depolarized  returns 
in  the  backscattering  direction  whereas  the  results  obtained  using  the  EBC  method  give 
significant  depolarized  returns  when  the  incident  direction  is  not  perpendicular  to  the  row 
direction  of  the  periodic  surface. 

We  have  also  derived  a  general  mixing  formula  for  discrete  scatterers  immersed  in  a 
host  medium.  The  inclusion  particles  are  assumed  to  be  ellipsoidal.  The  electric  field  inside 
the  scatterers  is  determined  by  quasistatic  analysis,  assuming  the  diameter  of  the  inclusion 
particles  to  be  much  smaller  than  the  wavelength.  The  results  are  applicable  to  general 
multiphase  mixtures,  and  the  scattering  ellipsoids  of  the  different  phases  can  have  different 
sizes  and  arbitrary  ellipticity  distribution  and  axis  orientation,  i.e.,  the  mixture  may  be 
isotropic  or  anisotropic.  The  resulting  mixing  formula  is  nonlinear  and  implicit  for  the 
effective  complex  dielectric  constant,  because  the  approach  in  calculating  the  internal  field 
of  scatterers  is  self-consistent.  Still,  the  form  is  especially  suitable  for  iterative  solution. 
The  formula  contains  a  quantity  called  the  apparent  permittivity,  and  with  different  choices 
of  this  quantity,  the  result  leads  to  the  generalized  Lorentz  -  Lorenz  formula,  the  generalized 
Polder  -  van  Santen  formula,  and  the  generalized  coherent  potential  -  quasicrystalline 
approximation  formula.  Finally,  the  results  are  applied  to  calculating  the  complex  effective 
permittivity  of  snow  and  sea  ice. 
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The  study  of  the  strong  fluctuation  theory  for  a  bounded  layer  of  random  discrete 
scatterers  are  extended  to  include  high-order  co-polarized  and  cross- polarized  second  mo¬ 
ments.  The  backseat tering  cross  sections  per  unit  area  are  calculated  by  including  the 
mutual  coherence  of  the  fields  due  to  the  coincidental  ray  paths  and  that  due  to  the  op¬ 
posite  ray  paths  which  are  corresponding  to  the  ladder  and  cross  terms  in  the  Feynman 
diagrammatic  representation.  It  is  proved  that  the  contributions  from  ladder  and  cross 
terms  for  co-polarized  backscattering  cross  sections  are  the  same,  while  the  contributions 
for  the  cross- polarized  ones  are  of  the  same  order.  The  bistatic  scattering  coefficients  in 
the  second-order  approximation  for  both  the  ladder  and  cross  terms  are  also  obtained.  The 
enhancement  in  the  backscattering  direction  can  be  attributed  to  the  contributions  from 
the  cross  terms. 

We  have  derived  the  dyadic  Green’s  function  for  a  two-layer  anisotropic  medium. 
The  Born  approximation  is  used  to  calculate  the  scattered  fields.  With  a  specified  cor¬ 
relation  function  for  the  randomness  of  the  dielectric  constant,  the  backscattering  cross 
sections  are  evaluated.  The  analytic  expressions  for  backscattering  coefficients  are  shown  to 
include  depolarization  effects  in  the  single-scattering  approximation.  It  is  also  shown  that 
the  backscattering  cross  section  (per  unit  area)  of  horizontal  polarization  can  be  greater 
than  that  of  vertical  polarization  even  in  the  case  of  half-space.  The  bistatic  scattering 
coefficients  are  first  calculated  and  then  integrated  over  the  upper  hemisphere  to  be  sub¬ 
tracted  from  unity,  in  order  to  obtain  the  emissivity.  The  principle  of  reciprocity  is  then 
invoked  to  calculate  the  brightness  temperatures.  It  is  shown  that  both  the  absorptive 
and  randomly  fluctuating  properties  of  the  anisotropic  medium  affect  the  behavior  of  the 
resulting  brightness  temperatures  both  in  theory  and  in  actual  controlled  field  measure¬ 
ments.  The  active  and  passive  results  are  favorably  matched  with  the  experimental  data 
obtained  from  the  first-year  and  the  multiyear  sea  ice  as  well  as  from  the  corn  stalks  with 
detailed  ground-truth  information. 
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Electromagnetic  wave  propagation  and  scattering  in  an  anisotropic  random  medium 
has  been  examined  with  Dyson  equation  for  the  mean  field  which  is  solved  by  bilocal  and 
nonlinear  approximations  and  with  Bethe-Salpeter  equation  for  the  correlation  of  field 
was  derived  and  solved  by  ladder  approximation.  The  effective  propagation  constants  are 
calculated  for  the  four  characteristic  waves  associated  with  the  coherent  vector  fields  prop¬ 
agating  in  an  anisotropic  random  medium  layer,  which  lire  the  ordinary  and  extraordinary 
waves  with  upward  and  downward  propagating  vectors.  The  z-component  of  the  effective 
propagation  constant  of  the  upward  propagating  wave  is  different  from  the  negative  of  that 
of  the  downward  propagating  wave,  not  only  for  the  extraordinary  wave  but  also  for  the 
ordinary  wave.  This  is  due  to  the  tilting  of  the  optic  axis  which  destroyes  the  azimuthal 
symmetry. 

Since  both  snow  and  ice  exhibit  volume  scattering  effects,  we  model  the  snow-covered 
ice  fields  by  a  three-layer  random  medium  model  with  an  isotropic  layer  to  simulate  snow, 
an  anisotropic  layer  to  simulate  ice,  and  the  bottom  one  being  ground  or  water.  The  snow 
and  ice  are  characterized  by  different  dielectric  constants  and  correlation  functions.  The 
theoretical  results  are  illustrated  for  thick  first-year  sea  ice  covered  by  dry  snow  at  Point 
Barrow  and  for  artificial  thin  first-year  sea  ice  covered  by  wet  snow  at  CRREL.  The  radar 
backscattering  cross  sections  are  seen  to  increase  with  snow  cover  for  snow-covered  sea 
ice,  because  snow  gives  more  scattering  than  ice.  The  results  are  also  used  to  interpret 
experimental  data  obtained  from  field  measurements. 
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ABSTRACT 

The  two-layer  anisotropic  random  medium  model  is  used  to  investigate  fully  polari- 
metric  scattering  properties  of  earth  terrain  media.  The  polarization  covariance  matrices 
for  the  untilted  and  tilted  uniaxial  random  medium  are  evaluated  using  the  strong  fluc¬ 
tuation  theory  and  distorted  Bom  approximation.  In  order  to  account  for  the  azimuthal 
randomness  in  the  growth  direction  of  leaves  in  tree  and  grass  fields,  an  averaging  scheme 
over  the  azimuthal  direction  is  also  applied.  It  is  found  that  characteristics  of  terrain 
clutter  can  be  identified  through  the  analysis  of  each  element  of  the  covariance  matrix. 
Theoretical  results  are  illustrated  by  the  comparison  with  experimental  data  provided  by 
MIT  Lincoln  Laboratory  for  tree  and  grass  fields. 


I.  INTRODUCTION 


The  layered  random  medium  model  for  terrain  cover,  in  conjunction  with  the  ap¬ 
plication  of  electromagnetic  wave  theory,  provides  a  systematic  approach  in  relating  the 
radar  backscatter  response  to  physical  properties  of  geophysical  media,  where  the  volume 
scattering  and  anisotropic  effects  are  attributed  to  the  embedded  inhomogeneities  with 
elongated  geometry  and  preferred  alignment.  For  active  and  passive  microwave  remote 
sensing,  the  effectiveness  of  this  model  has  been  demonstrated  by  extraction  of  physical 
parameters  from  data  matching  for  scene  interpretation  and  feature  identification  [Zuniga 
et  al.,  1980;  Tsang  and  Kong,  1981;  Tsang  et  al.,  1982;  Lee  and  Kong,  1985;  Lin  et  al., 
1988a,  b].  In  polarimetric  microwave  remote  sensing,  the  random  medium  model  has  also 
been  applied  to  radar  image  simulation,  terrain  clutter  identification  and  classification, 
and  radar  range  profile  simulation  of  terrain  clutter  [Shin  et  al.,  1986;  Borgeaud  et  al., 
1987a;  Borgeaud,  1987b;  Kong  et  al.,  1988;  Yueh  et  al.,  1988a;  Swarts  et  al.,  1988]. 

The  theoretical  modeling  of  earth  terrain  media  is  described  in  Section  II.  To  account 
for  the  large  permittivity  contrast  between  the  background  medium  and  the  embedded  in¬ 
homogeneities,  the  strong  fluctuation  theory  is  applied  in  Section  III,  in  order  to  derive 
the  external  field  which  is  decomposed  into  mean  (coherent  part)  and  scattered  field  (in¬ 
coherent  part)  components.  The  mean  field  is  derived  using  the  Feynman  diagrammatic 
technique  and  bilocal  approximation.  In  Section  IV,  the  effective  permittivity  tensor  of 
the  medium  is  obtained  from  the  constitutive  relation.  The  effective  permittivity  tensor 
is  computed  in  the  low  frequency  limit  for  an  unbounded  uniaxial  random  medium  with 
two-phase  mixture.  The  distorted  Born  approximation  is  then  used  to  calculate  the  scat¬ 
tered  field  in  Section  V.  The  polarization  covariance  matrices  for  the  untilted  and  tilted 
uniaxial  random  medium,  as  well  as  that  obtained  from  the  azimuthal  averaging  of  the 
tilted  case  are  discussed  in  Section  VI.  Theoretical  calculations  of  the  fully  polarimetric 
random  medium  model  are  compared  with  the  experimental  data  for  trees  and  a  grass  field 
supplied  by  MIT  Lincoln  Laboratory. 

II.  TWO-LAYER  ANISOTROPIC  RANDOM  MEDIUM  MODEL 

The  two-layer  anisotropic  random  medium  model  is  applied  to  simulate  earth  terrain 
media  such  as  tree  canopy  and  a  grass  field.  In  region  1  of  Figure  1,  the  random  medium 
with  height  d  is  characterized  by  a  spatially  random  permittivity  tensor  fi(F).  Regions 
0  and  2  (air  and  ground)  are  considered  to  be  homogeneous  media  with  permittivities  eo 
and  «j,  respectively.  All  three  regions  are  assumed  to  have  the  same  permeability,  p 0- 
The  polarimetric  radar  system  is  located  in  region  0.  The  coordinate  system  (x,  y,  z)  is 


oriented  so  that  the  zy  plane  coincides  with  the  air/clutter  interface.  Interfaces  at  z  —  0 
and  z  =  —d  are  assumed  to  be  planar,  extending  to  infinity  on  the  horizontal  plane,  and 
parallel  to  each  other.  The  volume  scattering  effect  is  caused  by  randomly  distributed 
scatterers  embedded  in  the  layer  (e.g.,  moisture  content  in  leaves)  while  the  anisotropic 
effect  is  due  to  the  elongated  shape  and  preferred  alignment  of  the  scatterers.  As  shown  in 
Figure  2,  the  moisture  content  in  trees  and  grass  is  modeled  as  cylindrical  scatterers  with 
the  preferred  alignment  in  the  yz  plane.  The  elongated  direction  is  along  the  z'-axis  and 
tilted  by  an  angle  with  respect  to  the  z-axis.  Hence,  the  coordinate  system  (*',  y\  z') 
which  describes  the  orientation  of  the  scatterers  is  rotated  clockwise  by  an  angle  rj)  about 
the  x-axis.  In  this  manner,  the  z'-axis  can  be  chosen  as  the  optic  axis  of  the  permittivity 
tensor  t|(F)  [Kong,  1986]. 

In  the  polsurimetric  microwave  remote  sensing  of  terrain  cover,  the  polarization 
covariance  matrix  is  essential  for  terrain- feature  identification  and  classification  [Borgeaud 
et  al.,  1987a;  Kong  et  ai,  1988].  For  a  plane  wave  impinging  on  the  random  medium  layer, 
the  scattered  electric  field  E'0(f)  is  related  to  the  incident  electric  field  UW(F)  by 
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where  the  horizontal  and  vertical  components  of  the  the  incident  and  scattered  electric 
fields  are  expressed  as  Eht,  Eh,,  and  E,„  respectively.  For  reciprocal  media,  the 

relation,  A*  =  A*>  is  obtained  in  the  backscattering  direction.  Thus,  the  covariance 
matrix  is  defined  as 

(l/»l’>  (mu)  (mu)' 

(mu)  (\M)  (mu)  (2) 

(mu)  (mu)  (i/..r) . 

where  A  is  the  illuminated  area  and  ( X ^  denotes  the  ensemble  average  of  the  random 
variable  X.  The  covariance  matrix  is  normalized  so  that  the  diagonal  terms  represent 
the  backscattering  cross  sections  per  unit  area  for  the  HH,  HV,  and  VV  polarizations, 
respectively. 


HI.  STRONG  FLUCTUATION  THEORY 


Consider  a  linearly  polarized  time-harmonic  electromagnetic  plane  wave  impinging 
on  the  random  medium.  The  time-harmonic  factor  e"*"*  is  used.  The  propagation  vector 
Km  at  the  observation  point  has  the  incident  angle  with  respect  to  the  z-axis  (Fig. 
1)  and  the  azimuthal  angle  <f>m  .  The  total  electric  fields  in  regions  0  and  1  satisfy  the 
following  vector  wave  equations: 


V  x  V  x  ~E0{f)  -  k£E0(f)  =  0 

(3) 

V  x  V  x  E,(F)  -  k’J^  ■  E,(F)  =  0 

*0 

(4) 

where  k £  =  In  the  microwave  frequency  range,  the  dielectric  constants  of  air  and 

moisture  content  are  very  distinct.  In  order  to  account  for  the  large  permittivity  contrast 
between  the  host  medium  and  embedded  inhomogeneities,  a  deterministic  quantity  is 
introduced  in  the  strong  fluctuation  theory  [Ryzhov  and  Tamoikin,  1970;  Tsang  and  Kong, 
1981;  Tsang  et  al.,  1982;  Lin  et  al.,  1988a,  b].  Equation  (4)  can  then  be  rewritten  as 

V  X  v  X  E,(f)  -  kl‘-2-  ■  E,(f)  =  kl  (I,(V  ~  fo)  ■  E,(f) 

=  k$t(f)  ■  E,(F)  (5) 

Physically,  lfi  is  the  effective  permittivity  tensor  of  the  medium  in  the  low  frequency 
range  where  the  scattering  loss  is  negligible  because  the  size  of  scatterers  is  much  smaller 
than  the  wavelength  of  the  incident  field  [Tsang  et  al.,  1985].  Treating  the  term  on  the 
right-hand  side  of  (5)  as  the  effective  source,  the  total  electric  fields  in  regions  0  and  1  can 
be  represented  in  integral  forms  as  follows: 

Z0(F)  =  E* \f)  +  kl  /  •  Qx(rx)  •  (6) 

E,(F)  =  ET(F)  +  Jv  d,F15„(F,Fl)  •  5,(F,)  ■  E,(F,)  (7) 

where  the  unperturbed  electric  Held.,  E^*(F)  end  E|  (F),  ere  solution,  to  the  homoge- 
neous  vector  wave  equations  in  (3)  and  (5)  in  the  absence  of  the  effective  source  term  and 
the  subscript  Vi  denotes  the  volume  of  region  1 .  The  dyadic  Green’s  functions,  (7oi(F,  Ft) 
and  Un(?,  Fj),  which  correspond  to  the  responses  at  F  in  regions  0  and  1,  respectively, 
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due  to  a  point  so orce  at  rt  in  the  anisotropic  homogeneous  medium  with  the  permittivity 
tensor  !fl ,  are  governed  by  the  following  vector  wave  equations: 


V  x  V  x  <?oi(r,  t*i  )  -  Jb0J£oi(F,  rt)  =  0 


V  X  V  X  (?n(r,  rt)  -  •  SU(F,  Ft)  =  7$(F,  Fi) 


(8) 


(9) 


When  the  source  and  the  observation  points  are  in  the  same  region,  there  a  singularity 
in  the  dyadic  Green’s  function,  which  can  be  dealt  with  by  separating  (7U(7S  Ft)  into 
principal  value  PS&n(r,f t)  and  singular  components,  namely, 


3i 


Zrnto  ri)  =  PSX?u(r,  fi)  -  tj^(F  -  Ft) 


(10) 


Thus,  substituting  (10)  into  (7)  yields 

J\W  =  ElV)  +  *J  /  ■£,(?,)  ■7.(1'.)  (11) 

"V\ 

^i(?)  =  P  +  **?t(?)]  *^i(?)  (12) 

?i00  =  *  P  +  •  ^(F)]  1  (13) 

Physically,  Fx(r)  can  be  interpreted  as  the  external  field  [Ryzhov  and  Tamoikin,  1970]. 

The  external  field  can  be  decomposed  into  a  coherent  part  (mean  field) 
and  an  incoherent  part  ?i(F)  (scattered  field)  [Tsang  et  ai,  1985].  Under  the  bilocal 
approximation,  the  mean  field  ^7\(F)^  cam  be  expressed  as 

(7\(F))  =  ^(F)  +  kl  fjy  d,Fld*FJP5^n(F,F1)  -^(F^F,) '  (14) 


where 


(15) 


<W/(rl.rl)  =  *o({l(fl)'P®®u(rl'r<)'£l(rl)) 

Using  the  index  notation,  (15)  can  be  expressed  as 

pfu/rO'!.*)],  =  ^r^^Ft.F,)  pll(Fi,F,)]mw  +  Tlmn,smn6(rx  -  f,)  (16) 
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where  Equation  (10)  wu  used  and  ^£ih»(Fi)£i»*(Fj)^  has  been  replaced  with  the  prod* 

uct  of  the  variance,  rJtw.f  (f,m,n,p  =  *,y,  z),  and  the  normalised  correlation  function, 
Rf(¥ i,Fj).  For  a  statistically  homogeneous  medium,  the  normalized  correlation  function 
is  a  function  of  the  displacement  between  Ft  and  F,,  i.e., 

*c0ri.*s)*  Ri(rx-rt)  (17) 

Notice  that  the  intrinsic  properties  of  the  embedded  inhomogeneities  are  directly  related 
to  the  volume  scattering  mechanism  through  the  normalized  correlation  function  in  (14). 

IV.  EFFECTIVE  PERMITTIVITY 

The  effective  permittivity  tensor  for  the  uniaxial  random  medium  is  derived  in  this 
section,  following  the  usual  approach  where  the  effect  of  boundary  layers  is  neglected  [Tsang 
and  Kong,  1981;  Stogryn,  1984].  In  an  unbounded  uniaxial  medium,  with  a  quasi-static 
effective  permittivity  tensor  of  the  form 

et  0  0  ' 

!fl=  0  «,  0  (18) 

.  0  0  tj, . 

the  mean  field  satisfies  the  following  vector  wave  equation 

V  x  V  x  {?»)  -kff-  (7,(F))  +  kl  <J*F,(„„(F, F.)  •  (l\(F.))  =  0  (19) 

The  dispersion  relation,  obtained  from  (19),  allows  the  effective  permittivity  tensor  for  the 
medium  to  be  defined  as 

«i  0  0  ■ 

fffi)  -  0  ti  0  (20) 

.0  0  <u ■ 

where  £Uff(k)  is  the  Fourier  transform  of  —  Fj)  as  given  in  (16).  In  the  low 

frequency  limit,  can  be  approximated  by 

[{, .„(»)]„  =  /’  <«*«(*)  p.(*)]_  +  (2i) 
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where  $<(£)  is  the  Fourier  transform  of  the  normalized  correlation  function  f?((r).  The 
Fourier  transform  of  the  dyadic  Green’s  function  Vt{f)  is  given  as 

r  o 

0 

0. 


*.(J) = 


(*; +*:-  *j«,)  k 


-k,km  kl 

0  0 


*; +?■(*!-  *;«..) 


r  (*;  -  »,.)  *.*,  (*;  -  *.■«,.) 


Jfc» 


*,*.  (*;  -  (*;  -  *:«,.) 


Jb> 


jfcj 


k.k. 


k.k, 

k,k. 

kl  ~  K  e. 


(22) 

in  which  k*  =  At*  +  fc*.  Analytical  expressions  for  {w/(0)  and  the  diagonal  tensor  3?x 
are  shown  in  Appendix  A.  For  a  two-phase  mixture  consisting  of  a  background  medium 
and  scatterers  with  permittivities  e»  and  €„  respectively,  ffl  can  be  evaluated  from  the 


following  criterion  [Tsang  and  Kong,  1981]: 

o 

II 

£ 

(23) 

From  (13)  and  (23),  we  obtain 

*+V V-'->+wVV-° 

(24) 

(25) 

where  f,  is  the  fractional  volume  of  scatterers.  Equations  (24)  and  (25) 

are  the  anisotropic 

version  of  Polder  and  van  Santen’s  mixing  formula  [Polder  and  van  Santen,  1946]. 


V.  DISTORTED  BORN  APPROXIMATION 

In  the  strong  fluctuation  theory,  the  original  random  medium  is  replaced  with  an 
equivalent  random  medium  which  has  the  permittivity,  ffl.  This  can  be  accomplished 
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without  altering  the  characteristic  of  the  volume  scattering  effect  because  the  effective 
source  defined  in  (5)  remains  intact,  that  is, 

*{  5.(1-)  •  -«?,(»-) -V.w  (26) 

where  (12)  and  (13)  have  been  used.  From  (6),  the  total  electric  field  in  region  0  can  be 
rewritten  as 

Zo(f)  =  >(r)  +  fco  /  <^F|f?*i(F,  Fv)  *  fl(?l)  •  7x(Fi) 

•Vi 

=  siV) + ir0(f)  (27) 

where  ^(f)  is  the  scattered  electric  field  due  to  the  effective  source.  By  adopting  the 
concept  of  the  distorted  Born  approximation  in  quantum  mechanics  [Newton,  1966;  Schiff, 
1968],  the  scatterers  are  assumed  to  be  embedded  in  the  equivalent  medium  with  the  effec¬ 
tive  permittivity,  fw//(0)  [Lang,  1981],  so  that  the  mean  field,  ^7\(F)^,  is  used  to  ap¬ 
proximate  the  external  field,  7*i (F),  and  the  mean  dyadic  Green's  function,  ^(?oi(F,  F,)^, 
is  used  to  replace  <?oi(F,  Ft).  The  following  vector  wave  equations: 

V  x  V  x  (Tx(7)}  -  •  (7\(F))  =  0  (28) 

V  x  V  x  (<?l0(F,  Fx)^  -  •  (<7io(f,  Ft)^  =  0  (29) 

in  conjunction  with  the  symmetric  property  of  the  dyadic  Green’s  functions,  namely, 
^01  (F,  rt  ,  are  used  to  derive  the  mean  field  and  the  mean  dyadic 

Green’s  function  for  the  scattered  electric  field  which  can  be  written  as 

zr„(?)  =  k}Jv  F,))  {,(7.)  •  (r, (!>,))  (30) 

Physically,  the  scattered  electric  field,  under  the  distorted  Born  approximation,  corre¬ 
sponds  to  the  single  scattering  of  the  coherent  fields  [Tsang,  et  al.,  1985].  It  is  also  known 
as  the  first-order  multiple  scattering  [Ishimaru,  1978].  After  decomposing  2?£(F)  into  hor¬ 
izontal  and  vertical  components  and  making  use  of  (1)  and  (2),  the  nine  elements  of  the 
polarization  covariance  matrix  can  be  computed. 


VI.  DATA  MATCHING  AND  APPLICATIONS 


Data  Matching 


For  an  untilted  uniaxial  random  medium,  the  optic  axis  is  in  the  z-direction,  V’  = 
=  0.  For  vertically  and  horizontally  polarized  incident  fields,  the  single  scattering 
process  in  the  distorted  Bora  approximation  does  not  depolarize  the  incident  field  in  the 
backscattering  direction  and  the  coefficient  fk9  vanishes.  Hence,  the  covariance  matrix 
becomes 


i  o 
0  0  0 


(31) 


UV 7  0  7  J 

where  o  is  the  backscattering  coefficient  for  the  HH  polarization,  7  is  the  copolarization 
intensity  ratio  (cr,./crkA),  and  p  is  the  normalized  correlation  coefficient  between  HH  and 
VV  returns  given  by 


(a./:.) 


(316) 


In  general,  for  the  case  of  tilted  uniaxial  random  medium  ^  0),  depo¬ 

larization  effects  exist  even  in  the  single  scattering  process,  and  the  nine  elements  of  the 
covariance  matrix  are  all  nonzero,  i.e., 


r  1  pyfi 

/3',/e  e  ty/=f t 
p'yfi  Cy/7*  7 


(32) 


where 


(32a) 

(326) 


are  the  correlation  coefficients  between  the  HH  and  HV  and  HV  and  W  channels,  respec¬ 
tively.  However,  when  the  downward  propagation  vector  lies  in  the  yi  plane,  i.e., 
<j>ot  =  90*,  both  the  double  refraction  phenomenon  and  the  depolarization  effect  in  the 
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backscattering  direction  disappear  so  that  /h  =0  [Lee  and  Kong,  1985a].  Under  this 
condition,  the  covariance  matrix  also  retains  the  same  form  as  the  untilted  case,  as  shown 
in  (31).  It  should  be  noted  that  the  growth  directions  of  leaves  in  trees  and  grass  fields  are 
not  necessarily  aligned  in  one  direction.  In  order  to  account  for  the  random  orientation 
of  the  scatterers  in  the  azimuthal  direction,  an  azimuthal  averaging  scheme  [Borgeaud, 
1987b]  is  applied  to  derive  the  polarization  covariance  matrix.  It  is  found  that  although 
fhm  ^  0,  four  of  the  elements  in  the  covariance  matrix  reduce  to  zero  [Borgeaud  et  al., 
1987a;  Kong  et  al.,  1988].  Thus  the  polarization  covariance  matrix  becomes 

’  1  0  P\fi' 

?  =  <r  0  .  0  (33) 

■PV 7  0  7  • 

The  existence  of  the  four  zero  elements  is  due  to  the  fact  that  the  random  medium  possesses 
an  azimuthal  symmetry  in  a  statistical  sense.  The  four  zero  elements  of  the  covariance 
matrix  indicate  that  there  is  no  correlation  between  the  copolarisation  (HH  and  W)  and 
cross- polarization  (HV)  returns,  i.e.,  /?  =  0  and  (  =  0. 

In  order  to  illustrate  the  random  medium  model,  polarimetric  data  for  forests  and 
grass  fields  obtained  from  MIT  Lincoln  Laboratory  are  matched  with  this  model.  The  op¬ 
erating  frequency  was  35  GHz,  and  the  angle  of  incidence  was  82*  .  Measured  covariance 
matrices  for  grass  and  tree  regions,  given  in  Table  1,  were  calculated  from  the  blocked 
out  regions  shown  in  Figure  7  of  [Swartz  et  al,  1988].  Analysis  of  this  database  indi¬ 
cated  that  there  is  essentially  no  correlation  between  the  HH  and  HV,  and  between 
the  HV  and  VV  polarimetric  returns;  from  a  statistical  point  of  view  the  terrain  clutter 
exhibited  azimuthal  symmetry,  and  therefore  0^0  and  {  ^  0  .  Thus,  the  form  of  the 
covariance  matrix  shown  in  (33)  was  utilised.  The  corresponding  theoretically  calculated 
covariance  matrix  parameters,  also  shown  in  Table  1,  were  obtained  by  azimuthal  averag¬ 
ing  [Borgeaud,  1987b]  of  the  covariance  matrices  obtained  using  the  two-layer  anisotropic 
random  medium  model  with  the  application  of  the  strong  fluctuation  theory  and  distorted 
Born  approximation. 

Applications 

The  random  medium  model  provides  a  systematic  approach  for  the  design  of  optimal 
target  detection  and  classification  algorithms  [Novak  et  al.,  1987]  and  the  identification  of 
terrain  media  such  as  vegetation  canopy,  forest,  and  snow-covered  fields  using  the  optimum 
polarimetric  classifier  [Kong  et  al.,  1988].  The  polarization  covariance  matrices  for  various 


terrain  corera  were  computed  from  theoretical  models  of  random  medium.  The  optimal 
classification  scheme  made  use  of  a  quadratic  distance  measure  and  was  applied  to  classify 
a  vegetation  canopy  consisting  of  both  trees  and  grass;  i.e.,  using  the  fully  polarimetric 
covariance  matrix  for  the  scattered  fields,  the  Bayes  likelihood  ratio  test  was  performed  on 
specific  measurements  to  classify  the  terrain  into  different  categories.  The  Bayes  likelihood 
ratio  test  has  been  shown  to  be  optimal  in  the  sense  that  it  minimises  the  probability  of 
error  [Fukunaga,  1972;  Hord,  1986;  Richards,  1986;  Swain,  1978]. 

The  probabilities  of  error  obtained  using  fully  polarimetric  data  were  compared  with 
the  probabilities  of  error  obtained  using  several  single,  polarimetrically  derived  features. 
For  each  single  feature  studied,  the  corresponding  probability  density  function  was  derived. 
Once  the  probability  density  functions  were  known,  the  Bayes  likelihood  ratio  test  was 
performed  to  classify  terrain  elements  into  different  categories.  For  classification  schemes 
based  on  single  features,  closed-form  expressions  for  the  probability  of  error  were  calculated 
[Kong  et  al.,  1988]. 

A  supervised  Bayes  classification  was  applied  to  synthetic  aperture  radar  (SAR) 
polarimetric  images  in  order  to  identify  their  various  earth  terrain  components  [Lim  et 
al.,  1988].  Both  fully  polarimetric  and  normalised  polarimetric  classifications  were  em¬ 
ployed  to  classify  radar  imagery.  It  was  again  shown,  in  this  case  through  use  of  radar 
images,  that  fully  polarimetric  classifications  yielded  optimal  results;  however,  an  optimal 
normalised  classification  scheme  [Yueh  et  al.,  1988b]  indicated  improved  performance  in 
regions  where  the  absolute  backscattering  coefficients  differed  significantly  from  that  of  the 
training  regions  due  to  the  variation  in  return  power  as  a  function  of  incident  angle. 

The  Bayes  classification  is  known  to  give  minimum  probability  of  error  as  long  as 
the  statistical  distributions  of  the  radar  returns  are  known  and  the  training  areas  are 
accurate.  However,  as  noted  above,  the  selected  training  regions  were  not  sufficient  to 
classify  correctly  the  entire  image  due  to  the  variation  in  incident  angle  over  the  imaged 
swath.  By  employing  the  random  medium  model  to  generate  the  classifier  training  data, 
the  effect  of  the  incident  angle  can  be  included  in  the  classification  scheme. 

Another  application  of  the  fully  polarimetric,  multi-incident  angle  and  frequency 
random  medium  model  is  to  match  training  data  from  various  image  regions,  say  at  one 
frequency  and  during  one  particular  season,  and  then  generate  classifier  training  data 
which  will  predict  the  backscattered  response  of  the  terrain  at  different  frequencies  and 
during  different  seasons  of  the  year.  For  example,  the  effects  of  the  winter  season  can  be 
simulated  by  adding  a  layer  of  snow  to  a  terrain  from  which  training  data  was  recorded 


during  the  summer  season.  In  this  manner,  data  need  not  be  collected  for  all  seasons  and 
at  all  frequencies.  The  random  medium  model  can  be  used  to  predict  the  backscattered 
response  based  on  the  initial  training  data. 


APPENDIX  A 

The  normalised  correlation  function,  y\  z'),  for  the  unbounded  uniaxial 

random  medium  in  the  (i1,  y*,  z')  coordinate  frame  is  assumed  to  be 


exp 


K  K  is  J 


so  that  is  given  by 


where 


and 


T  ,„.(/  + S)  0  0 

fw,(0)  =  0  ^(1  +  5)  0 

0  0  r  „„(/,  + s.) 

r.«««  =  rtjrff  = 


r  _  *><  \  2  ft^+v2rn)(OT+i) 

^(a+2)lv^+Tog[  v^-iv/c^TTj 


„  /^C,  IT  1 

u+2log[  vm +••>/?  J  2} 


f  • gg  -  X) , 

™.(0  +  2)  \  iog 


{y/P^X  +  VFTl)  (#TI  +  l)' 

y/$-X-iyJx{P  +  l) 


(Al) 


(A2) 

(A3) 

(A4) 


+  2^(0  -  x)(x  +  2)  log 


V5(v^+2-l) 

,  *(x  +  2)l 

Vx  +  2  +  Vx  _ 

+  2  I 

(A5) 


where 


2<o  f  2(3  [(v^x  +  V^+T)(>^fT+l) 

**»•(&  +  2)  \v^TT  °®[  V/T^-iy/^TT) 

4  J [EE*  log  ^(VxT?-l)|  90  1 

V  x  +  2  °®  v'xTTJ  +  Vx  ~TJ 
<->«.  X  =  ^C,  a  =  <  +  5  *,&•« 


The  components  of  ^  are  obtained  from  the  frequency-independent  part  of  (A5)  and 
(A6): 


g - «<>  f  2rl*  ,  _  OvV  +  1  +  qJ 

^i  +  2)UJtfJT)lo*[  V^TT-l 


-4S&M 

5,  = ^ ]rir  bvV +  !  +  ?» 

e»*(^*  +  2)  \  n  J(„»  +  vV  +  1-1 


2%/2q  f>/2  +  ll  ,1 

-  °*  Ivf^rJ  ~ , } 


n  =  ^ 

’  >  V? 


(A10) 


From  (A8)  and  (A9),  it  can  be  shown  that  5  and  5,  satisfy  the  condition  [Stogryn 
1983): 

2is+‘ts-  = 1  <"« 
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Figure  1  :  Anisotropic  two-layer  random  medium  configuration. 


z 


Figure  2  :  Geometry  of  the  permittivity  tensor  in  the  anisotropic  random  medium. 
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Table  1  :  Covariance  matrix  elements  for  grass  (a)  and  trees  (b). 
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Abstract 


Fully  polarimetric  scattering  of  electromagnetic  waves  from  earth  terrain  media  are 
studied  with  a  three-layer  configuration.  The  inhomogeneous  layers  are  modeled  as  ran¬ 
dom  media  and  characterized  with  three-dimensional  correlation  functions.  The  strong 
fluctuation  theory  and  the  distorted  Bora  approximation  are  used  to  calculate  the  covari¬ 
ance  matrix  which  describes  the  polarimetric  radar  scattering  from  the  remotely  sensed 
media.  Theoretical  results  are  first  obtained  for  two-layer  cases  and  then  used  to  compare 
with  results  for  three-layer  cases  to  identify  the  effects  on  polarimetric  wave  scattering 
due  to  the  top  layer.  A  covering  dry-snow  layer  can  give  rise  to  the  enhancement  of  radar 
returns  from  first-year  sea  ice  and  also  the  oscillation  in  the  real  and  imaginary  parts  of  the 
correlation  coefficient  p  between  the  W  and  HH  returns  as  a  function  of  incident  angle. 
The  magnitude  of  p,  however,  does  not  show  the  oscillation  and  clearly  retains  the  similar 
characteristics  as  observed  directly  from  bare  sea  ice.  When  the  covering  layer  is  thick  or 
lossy,  some  polarimetric  information  from  the  sea-ice  layer  is  masked  out. 


I.  Introduction 


Polarimetric  radar  scattering  from  bare  terrain  fields  such  as  snow,  ice,  and  veg¬ 
etation  has  been  successfully  modeled  with  the  two-layer  random  medium  configuration 
[1-3].  In  this  paper,  the  random  medium  model  is  extended  to  a  three-layer  configuration 
to  account  for  fully  polarimetric  scattering  from  geophysical  media  under  the  effects  of 
weather,  diurnal  and  seasonal  variations,  and  atmospheric  conditions  such  as  ice  under 
snow,  meadow  under  fog,  and  forest  under  mist.  The  top  scattering  layer  is  modeled  as 
an  isotropic  random  medium  which  is  characterized  by  a  scalar  permittivity.  The  middle 
scattering  layer  is  modeled  as  an  anisotropic  random  medium  with  a  symmetric  permittiv¬ 
ity  tensor  whose  optic  sods  can  be  tilted  due  to  the  preferred  alignment  of  the  embedded 
scatterers.  The  bottom  layer  is  considered  as  a  homogeneous  half-space.  Volume  scattering 
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effects  of  both  random  media  are  described  by  three-dimensional  correlation  functions  with 
variances  and  correlation  lengths  corresponding  to  the  strengths  of  the  permittivity  fluctu¬ 
ations  and  the  physical  sizes  of  the  inhomogeneities,  respectively.  The  strong  fluctuation 
theory  [4]  is  used  to  derive  the  mean  fields  in  the  random  media  under  the  bilocal  approx¬ 
imation  with  singularities  of  the  dyadic  Green’s  functions  properly  taken  into  account  and 
effective  permittivities  of  the  random  media  are  calculated  with  two-phase  mixing  formu¬ 
las.  The  distorted  Bom  approximation  is  then  applied  to  obtain  the  covariance  matrix 
which  describes  the  fully  polarimetric  scattering  properties  of  the  remotely  sensed  media. 
The  polarimetric  information  proves  to  be  useful  in  the  identification,  classification,  and 
radar  image  simulation  of  earth  terrain  media. 

To  observe  polarimetric  scattering  directly  from  terrain  media  such  as  snow,  ice,  and 
vegetation,  the  three-layer  configuration  is  first  reduced  to  two-layers  by  removing  the  top 
layer.  Such  media  exhibit  reciprocity  as  experimentally  manifested  in  the  close  proximity  of 
the  measured  backscattering  radar  cross  sections  <7,*  and  <rk,  and  theoretically  established 
in  the  random  medium  model  with  symmetric  permittivity  tensors.  The  theory  is  used 
to  investigate  the  signatures  of  isotropic  and  anisotropic  random  media  on  the  complex 
correlation  coefficient  p  which  is  helpful  in  terrain  classification  [5].  Consideration  will  also 
be  given  to  the  nonzero  depolarization  covariance  terms  which  are  attributive  to  the  tilt 
of  the  optic  axis  in  the  anisotropic  scattering  layer. 

The  effects  on  polarimetric  wave  scattering  due  to  the  top  layer  are  identified  by 
comparing  the  three-layer  results  with  those  obtained  from  the  two-layer  configuration. 
The  theory  is  used  to  explain  the  enhancement  of  the  radar  returns  due  to  a  dry-snow 
cover  on  top  of  first-year  sea  ice,  predict  the  effect  of  the  snow  layer  on  p,  and  study  the 
masking  effects  on  the  radar  scattering  from  the  ice  layer.  In  this  application  of  the  random 
medium  model  to  polarimetric  remote  sensing  of  earth  terrain,  the  encountered  media  are 
reciprocal  and  characterized  by  symmetric  permittivity  tensors.  In  other  applications 
involving  a  nonreciprocal  medium,  the  permittivity  tensor  will  not  be  symmetric  as  in  the 
case  of  the  irregular  ionosphere  which  has  been  modeled  as  a  gyrotropic  random  medium 

[Si¬ 


ll.  Formulation 


The  scattering  configuration  is  depicted  in  Figure  1.  An  electromagnetic  wave 
is  incident  from  air  (region  0)  to  the  layered  media.  Region  1  is  a  scattering  medium 
with  isotropic  scatterers  randomly  embedded  such  as  snow,  fog,  or  mist  whose  electrical 
property  can  be  characterized  with  a  spatially  dependent  permittivity  Region  2  is 

an  anisotropic  random  medium  containing  nonspherical  scatterers  such  as  sea  ice,  grass,  or 
trees  which  has  a  spatially  dependent  permittivity  e»(F).  Region  3,  which  can  be  sea  water 
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or  ground,  is  considered  homogenous  with  a  deterministic  permittivity  e«.  The  infinite 
planar  interfaces,  the  thickness  of  the  layers,  and  the  coordinate  system  (z,  y,  z)  are  shown 
in  Figure  1.  Due  to  the  preferred  alignment  of  the  nonspherical  scatterers,  region  2  is 
effectively  uniaxial  with  the  optic  axis  z'  tilted  off  the  z-axis  by  an  angle  ^  in  the  yz-plane 
as  illustrated  in  Figure  2. 

The  horizontal  and  vertical  components  (Eki  and  Eti,  respectively)  of  the  incident 
electric  field  can  be  related  to  the  horizontal  and  vertical  components  ( Ek ,  and  2?„,  re¬ 
spectively)  of  the  scattered  electric  field  by  the  scattering  matrix  in  the  following  manner 
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In  the  polarimetric  remote  sensing  of  earth  terrain,  the  encountered  media  are  reciprocal. 
The  reciprocity  renders  /.*  =  fk, .  In  this  case,  the  fully  polarimetric  radar  backscattering 
from  the  layered  random  medium  is  characterized  by  a  three  by  three  covariance  matrix 
defined  by 
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where  A  is  the  illuminated  area,  the  asterisks  denote  the  complex  conjugates,  and  the 
angular  brackets  represent  the  ensemble  averages.  In  terms  of  the  polarimetric  radar 
backscattering  coefficients,  the  covariance  matrix  (2)  can  be  written  as 
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Note  that  the  diagonal  terms  *,»**,  and  <r„„  are  the  conventional  radar  backscat¬ 

tering  coefficients  <rHht  <r,fc,  and  <r„,  respectively.  Normalized  to  <r  =  cr**,  the  covariance 
matrix  can  be  formalized  as 
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where  e  and  7  are  the  intensity  ratios 
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and  the  correlation  coefficients  p,  0,  and  £  are 
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From  the  above  definitions,  the  polarimetric  radar  backscattering  coefficients  are 
determined  by  the  incident  and  scattered  electric  fields  as  follows  [7] 


4*r* 
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where  p,  u,  r,  or  k  can  be  h  or  v,  and  the  subscripts  i  and  s  denote  the  incident  and 
scattered  fields,  respectively.  According  to  (7),  the  covariance  matrix  can  be  calculated 
when  the  scattered  field  is  obtained. 


III.  The  Scattered  Field 


Consider  a  linearly  polarized  electromagnetic  plane  wave  incident  on  the  three-layer 
random  medium.  For  time-harmonic  waves,  the  total  fields  1?0(F),  Ia(F),  and  U,(F)  in  the 
respective  regions  0,  1,  and  2  satisfy  the  following  wave  equations 

V  x  V  x  Eo(F)  -  fco^o(F)  =  0  (8) 

V  x  V  x  ^(F)  -  =  0  (9) 

V  x  V  x  E,(F)  -  fc*  ^  •  E,(F)  =  0  (10) 

In  the  remote  sensing  of  earth  terrain,  a  strong  permittivity  contrast  between  the  scat- 
terers  and  the  background  medium  is  often  encountered.  The  strong  fluctuation  theory 
is  therefore  necessitated  in  the  calculation  of  the  random-medium  effective  permittivity. 
Deterministic  permittivities  egl  and  f,j  are  respectively  introduced  in  both  sides  of  (9)  and 
(10)  and  the  following  vectors 
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(13) 


are  treated  as  the  effective  sources  so  that  the  wave  equations  (9)  and  (10)  become 

vxvx  r,(F>  -  =  *0>5,(f)  •  Ei(f) 

*0 

V  X  V  X  E,(f)  -  *0’^  •  E,(f)  =  •  E,(f)  (14) 

The  mean  fields  in  the  random  media  are  then  derived  under  the  bilocal  approximation. 
In  the  bilocal  approximated  Dyson’s  equations,  the  possible  coincidence  of  an  observation 
point  with  a  source  point  gives  rise  to  singularities  in  the  dyadic  Green’s  functions  (DGFs). 
The  singularity  problem  is  resolved  by  decomposing  the  DGFs  into  principal  value_parts 
with  appropriate  exdusionjrolumes^  and  Dirac  delta  parts  with  dyadic  coefficients  and 
?).  The  quantities  <fl,  ffJ,  and  are  obtained  with  two-phase  mixing  formulas  and  the 
requirement  of  secular- term  elimination  [4].  The  variances  of  the  permittivity  fluctuations 
and  the  effective  permittivities  and  I w  for  regions  1  and  2  are  calculated  with  the 
Fourier  transform  method  for  a  given  correlation  function  [3,4]. 

With  the  obtained  effective  permittivities  which  account  for  both  the  dissipation 
and  the  scattering  losses,  the  scattered  field  under  the  distorted  Bora  approximation  is 

E.(F)  =  *o  fVt  **(&.(»',  ?i))  •  ?»(*)  •  (7 i(*i)) 

+  *0  fYt  f,))  •  *,(Ft)  •  (y*(F,))  (15) 

and  the  scattered  intensity  is 

(e.(f)  •  £>))  =  E  ko  /  <fr»  /  ^  -  *») 

•  ((c?oi«(F,Ft))(Fu(Ft)))  •  ((<?  Ol»l(F,  *>)(*-<*>»*  ■ 

+  E  kQ  f  C’tJiWm  (F|  ~  Fj) 

•  ((<WF, *))(«.(*)))  •  ((c«1(F,F,))(l«l',)))*  (16) 

where  Vi  and  Vj  are  the  spaces  respectively  occupied  by  regions  1  and  2  and  the  mean  DGFs 
(<?oi(F,Fn))  and  (Rtf,*.))  for  n  =  1,2  are  calculated  with  the  effective  permittivities. 

In  (15)  and  (16),  the  external  fields  7\(F)  and  7,(F),  the  scatterers  £t(F)  and  ?,(F),  and 
the  correlation  functions  C( i>»im(F)  and  C(j,«m(F)  are  defined  as 

^.(F)  =  [7  + 1,  •  5»(F)]  •  rn(F)  ,  n  =  1,2  (17) 
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U^)  =  $.(?)  •  p  +  ?.  *  $.(*)]  "l ,  n  =  1, 2  (18) 

<?<-,««  fa  ~  *)  =  (&;*fa)  &»fa))  .  n  =  1,2  (19) 

In  this  paper,  the  following  forms  of  the  normalized  correlations  are  used 

For  region  1  :  /2ei(r)  =  exp  _  M  — ^  (20) 

For  region  2  :  R(j(f)  =  exp  exp  (21) 

The  results  above  are  applied  in  (7)  to  obtain  the  covariance  matrix  of  the  forms  (3)  and 

(4). 

IV.  Results  and  Discussions 


The  three-layer  configuration  can  be  reduced  to  two-layers  by  removing  the  top  layer 
to  investigate  the  polarimetric  backscattering  directly  from  uncovered  earth  terrain  media. 
Consider  first  an  electromagnetic  wave  of  9  GHz  incident  on  an  isotropic  random  medium 
composed  of  a  background  permittivity  e»  =  (3.15  +  *0.002)4,  and  a  4.5%-volume  fraction 
of  isotropic  scatterers  having  a  permittivity  eit  =  (38  +  *41)4,  with  a  correlation  length  of 
1  mm.  The  strong  fluctuation  theory  gives  a  variance  of  0.32  and  an  effective  permittivity 
=  (3.6  +  t0.093)e<,.  Underlying  the  random  medium  of  thickness  dj  =  1.65  m,  the 
bottom  layer  has  a  permittivity  ea  =  (45  +  *40)4,.  The  physical  parameters  required 
to  calculate  the  mean  backscattered  field  are  summarized  in  Figure  3(a).  In  the  case 
of  Figure  3(b),  the  random  medium  has  the  same  physical  parameters  as  in  the  case  of 
Figure  3(a)  except  that  it  is  effectively  anisotropic  with  correlation  lengths  tp  —  1  mm 
and  l,  =  3  mm  for  which  the  strong  fluctuation  theory  yields  a  variance  8,  =  0.22  and 
an  effective  permittivity  —  (3.52  +  *0.08)4,  across  the  optic  axis  and  8,  =  1.1  and 
tU}t,  =  (3.99  +  *0.458)4,  along  the  optic  axis.  For  both  the  isotropic  and  untilted  (V*  =  0) 
anisotropic  random  media,  the  depolarization  terms  £•*•«,  and  are  zero  under 

the  first-order  distorted  Born  approximation  rendering  the  covariance  matrix  of  the  form 
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In  all  of  the  following  results,  the  incident  wave  vector  fco  is  taken  to  be  parallel 
to  the  zz-plane  though  the  theory  can  be  used  for  any  arbitrary  incident  direction.  The 
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conventional  backscattering  coefficients  <7**  and  <r„  are  plotted  as  a  function  of  incident 
angle  in  Figures  4(a)  and  4(b)  for  the  two-layer  isotropic  and  untilted  anisotropic  cases, 
respectively.  As  observed  from  Figure  4,  <r„  crosses  over  *t  an  incident  angle  of  about 
25*  for  the  anisotropic  case  whereas,  for  the  isotropic  case,  is  higher  than  over  the 
range  of  incident  angles  under  consideration.  The  distinction  of  the  conventional  backscat- 
tering  coefficients  <r„  and  crKk  associated  with  the  two  different  random  media  is,  however, 
not  as  obvious  as  that  of  the  polarimetric  correlation  coefficients  p  as  shown  in  Figure  5 
where  the  isotropic  random  medium  (a)  results  in  p  with  the  value  of  approximately  1.0 
over  the  concerned  incident  angles  and  the  untilted  anisotropic  random  medium  (b)  re¬ 
sults  in  p  with  complex  values  decreased  as  the  incident  angle  is  increased.  The  correlation 
coefficient  p  also  contains  information  about  the  tilt  of  the  optic  axis  in  the  anisotropic 
random  medium  as  illustrated  in  Figure  6(a)  where  the  maximum  magnitude  of  p  is  at 
normal  incident  in  the  untilted  case  (V»  =  0*)  and  moves  to  a  larger  incident  angle  as  the 
tilted  angle  becomes  larger.  Also  shown  in  Figure  6(b),  the  phase  of  p  does  not  change 
sign  for  =  0*  whereas,  in  the  tilted  cases,  the  phase  of  p  flips  sign  at  the  incident  angle 
where  the  magnitude  of  p  is  maximum. 

It  should  be  noted  that  the  tilt  of  the  optic  axis  is  also  related  to  the  nonsero 
depolarization  terms  in  the  covariance  matrix  which  reduce  to  zero  for  an  isotropic  or 
untilted  anisotropic  random  medium  under  the  Born  approximation  as  indicatedm  (22). 
For  the  purpose  of  comparison,  shown  below  are  the  covariance  matrices  (7.,  (7.10,  and 
calculated  at  an  incident  angle  of  40*  respectively  for  the  isotropic,  the  untilted,  the 
10*-tilted,  and  the  20*-tilted  anisotropic  random  media 
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To  identify  the  effects  of  the  covering  top  layer  on  the  radar  backscattering  from 
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the  lower  layer,  the  components  of  covariance  matrices  are  compared  for  the  two-layer  and 
the  three-layer  configurations.  In  Figure  7(a),  the  top  layer  is  a  low-loss  isotropic  random 
medium  of  thickness  dx  =  0.1  m  composed  of  a  background  of  permittivity  «u  =  eo  and 
a  25%-volume  fraction  of  isotropic  scatterers  having  a  permittivity  «i,  =  (3.15  +  *0.002)e<> 
with  a  correlation  length  it  =  0.3  mm.  From  the  strong  fluctuation  theory,  the  variance 
and  the  effective  permittivity  are  found  to  be  S  =  0.270  and  e^f  =  (1.38  +  *4.38  x  10'4)<o» 
respectively.  In  Figure  7(b),  the  top  layer  is  a  lossy  isotropic  random  medium  with  the 
same  physical  parameters  as  in  the  case  of  Figure  7(a)  except  that  the  imaginary  part  of 
the  scatterer  is  now  0.3  <o  which  results  in  S  =  0.275  and  t^}  —  (1.38  +  *0.039)eo  due  to 
the  strong  fluctuation  theory.  In  both  cases,  the  bottom  layer  and  the  middle  anisotropic 
random  medium  with  a  10#-tilted  angle  are  kept  the  same  as  in  the  case  of  Figure  3(b). 
Displayed  in  Figure  8  are  the  plots  of  <7^  and  <r„  as  a  function  of  incident  angle  for  the 
two-layer,  the  low-loss,  and  the  lossy  three-layer  configurations.  In  the  low-loss  case,  both 
<Tkh  and  <r„  are  enhanced.  Furthermore,  the  boundary  effect  is  manifested  in  the  form  of 
the  oscillation  on  <7**  and  <r„.  The  oscillation  can  also  be  seen  on  the  real  and  imaginary 
parts  of  the  correlation  coefficient  p  as  shown  in  Figure  9(a).  The  magnitude  of  p,  however, 
does  not  exhibit  the  oscillation  while  clearly  retaining  the  same  characteristics  as  observed 
directly  from  the  two-layer  configuration.  Thus,  the  correlation  coefficient  p  can  carry 
information  from  both  the  covering  low-loss  isotropic  layer  and  the  lower  tilted  anisotropic 
layer  in  a  rather  distinctive  manner.  In  contrast  to  the  low-loss  case,  the  lossy  top  layer 
can  diminish  both  <ru  and  <r„  and  depress  the  boundary-effect  oscillation  as  demonstrated 
in  Figure  8.  When  the  thickness  of  the  lossy  top  layer  increases,  the  behavior  of  the 
correlation  coefficient  p  becomes  more  and  more  similar  to  the  isotropic  case  (compare 
Figures  9(b)  with  5(a)).  In  that  sense,  the  lossy  top  layer  can  mask  out  information  from 
the  lower  anisotropic  layer. 

As  illustrated  in  this  section,  the  polarimetric  information  conveyed  in  the  covari¬ 
ance  matrix  is  important  in  the  identification  and  classification  of  earth  terrain  media. 
Actually,  the  parameters  chosen  in  the  two-layer  configuration  of  Figure  3(b)  physically 
characterize  a  sea-ice  layer  over  sea  water.  Due  to  the  physical  nature,  the  random  medium 
model  can  help  remotely  identify  the  ice  type  with  a  polarimetric  radar.  Also  in  Figure 
7(a),  the  physical  parameters  correspond  to  a  dry-snow  layer  covering  the  sea-ice  layer 
floating  on  sea  water.  In  this  case,  it  has  been  shown  above  that  the  polarimetric  volume 
scattering  effects  can  reveal  useful  information  from  the  lower  layer  even  if  it  is  covered 
under  another  layer.  This  is  helpful  in  the  remote  sensing  of  the  arctic  regions  in  the  clas¬ 
sification  of  snow-covered  ground  and  snow-covered  sea-ice  regions.  Since  the  grain  size 
of  the  snow  cover  is  modeled  with  the  correlation  length  [8],  the  theory  can  account  for 
diurnal  changes  which  affect  the  size  of  snow  grains.  In  the  remote  sensing  of  vegetation, 
the  three-layer  random  medium  model  is  very  versatile  because  it  can  include  the  effects 
of  weather  and  seasonal  changes  such  as  meadow  under  fog  or  forest  under  mist.  The 
examples  in  this  section  serve  the  illustration  of  the  physical  nature,  the  usefulness,  and 
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the  versatility  of  the  three-layer  random  medium  model  in  the  polarimetric  remote  sensing 
of  earth  terrain. 
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Figure  3  :  Physical  parameters  for  the  two-layer  isotropic  and  anisotropic  media. 
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Figure  4  :  Comparisons  of  the  conventional  backscattering  coefficients  <rhk  and  <r„  for 
the  two-layer  isotropic  and  untilted  anisotropic  cases. 
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igure  5  :  Comparison  of  the  correlation  coefficients  p  for  the  two-layer  isotropic  and 
anisotropic  cases.  In  (b),  the  left  vertical  scale  is  for  \p\  and  Rt{p)  and  the 
right  vertical  scale  is  for  Im{p}. 
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Figure  6  :  Magnitudes  and  phases  of  p  for  the  two-layer  anisotropic  cases. 
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Figure  7  :  Physical  parameters  for  the  three-layer  configurations  with  low-loss 
and  lossy  top  layers. 
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Figure  8  :  Comparisons  of  the  conventional  backscattering  coefficients  <rkk  and  <r„  for 
the  three-layer  10®-  tilted  anisotropic  cases  with  low-loss  and  lossy  top  layers. 
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Figure  9  :  Correlation  coefficients  p  for  the  three-layer  10®— tilted  anisotropic  cases  with 
low-loss  and  lossy  top  layers.  In  (a),  the  left  vertical  scale  is  for  |/?|  and 
Re{p}  and  the  right  vertical  scale  is  for  Im{p}. 


